Synthetic Mechanochemical Molecular Swimmer 
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A minimal design for a molecular swimmer is proposed that is a based on a mechanochemical 
propulsion mechanism. Conformational changes are induced by electrostatic actuation when specific 
parts of the molecule temporarily acquire net charges through catalyzed chemical reactions involving 
ionic components. The mechanochemical cycle is designed such that the resulting conformational 
changes would be sufficient for achieving low Reynolds number propulsion. The system is analyzed 
within the recently developed framework of stochastic swimmers to take account of the noisy en- 
vironment at the molecular scale. The swimming velocity of the device is found to depend on the 
concentration of the fuel molecule according to the Michaelis-Menten rule in enzymatic reactions. 
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One of the aims of nanotechnology is to be able to make 
synthetic molecular devices that could propel themselves 
through fluidic environments and perform targeted tasks 
such as deliveryof therapeutic agents or carrying out me- 
chanical work [ij. At such small scales, one cannot apply 
standard deterministic strategies used in engineering at 
the macroscopic scales, as the dynamics of any device 
will be overwhelmed by fluctuations of thermal or other 
origins. The right strategy is to find a way to impose 
a bias on these fluctuations such that they will average 
out to our desired behavior, as can be learned from the 
example of biological molecular motors [2| . 

In addition to enduring the noisy environment, molec- 
ular swimmers would also need to overcome the funda- 
mental difficulty posed by the governing hydrodynamic 
rules at low Reynolds number conditions: they should 
utilize a non- reciprocal sequence of deformations to break 
the time reversal symmetry in their periodic motions and 
achieve swimming [3|. This is a nontrivial task if one is 
to use a minimal set of degrees of freedom J4l4l l| , which 
will most likely be what can be afforded in a human- 
engineered device. 

There has been a significant recent development along 
this line with the advent of a number of experimentally 
made micron-scale prototypes of such low Reynolds num- 
ber swimmers. Using magnetic actuation of micron-sized 
linked magnetic beads 1^ Il3| or manipulation of par- 
ticles using optical tweezers [1J|, it was demonstrated 
that low Reynolds number propulsion can be achieved at 
the micro-scale via non-reciprocal motion. In these ex- 
periments, the actuation mechanisms were externally en- 
forced (using oscillating magnetic fields or optical traps), 
which means that strictly speaking the devices were not 
self-propelled swimmers. Moreover, it is not clear if such 
actuation mechanisms could be scaled down to be ap- 
plicable to molecular devices. This means that we are 
in need of new strategies to be able to program non- 
reciprocal molecular deformations. 

While other external triggers — such as laser pulses that 



could induce conformational changes — could in principle 
be pursued py] , it would be ideal to be able to couple the 
conformational changes to a local source of free energy 
via a chemical reaction, like motor proteins [l5|. This 
will allow the mechanism to potentially accommodate 
arbitrary degrees of complexity, and avoid interference 
between different parts of the system. To my knowledge, 
this scenario has not been considered so far in the litera- 
ture for microscopic swimmers, although there have been 
demonstrations of mechanical actuation of elastic gels 
by using oscillatory chemical reactions (at macroscopic 
scales) [l6| . A strategy based on catalysis of a chemical 
reaction has already been used for making self-propelled 
colloidal particles, which utilize non-equilibrium inter- 
facial (phoretic) interactions towards their net motion 
[iTI . lis! . While this mechanism of motility has shown to 
be extremely powerful, exploring other possible strate- 
gies would help create a diverse range of approaches that 
could be used in the future to optimize devices for spe- 
cific purposes and overcome strategy-specific hurdles for 
a particular end use. For example, a swimmer that uti- 
lizes interfacial interactions as its main motility mech- 
anism is potentially vulnerable to intervention by other 
interfacial forces when it is close to boundaries or other 
devices. Moreover, one wonders why nature has not cho- 
sen to take advantage of the "phoretic" strategy in any 
living system; one may speculate that it could be because 
it cannot be conveniently miniaturized and robustly func- 
tion in maximally crowded intracellular environments. 

Here I propose an alternative minimal design for a 
molecular swimmer, which would be able to catalyze 
chemical reactions and use the free energy gain and prod- 
ucts to induce the conformational changes that would be 
sufficient for achieving low Reynolds number propulsion 
in a noisy environment [19| . The present study could also 
shed some light on how the motor proteins could com- 
bine enzymatic action with conformational changes, as it 
presents a full physical analysis of a simple model that 
explains how such a coupling could work in practice. 
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FIG. 1. The 6-state chemical cycle of the three-sphere molec- 
ular device that comprises enzyme h at the head and enzyme 
t at the tail (see Fig. [2)|. h catalyzes F — > Q~ + G^ and t 
catalyzes R~ -|- P"*" — > S. 



The proposed design for the molecular swimmer is 
based on the sirnple three-sphere low Reynolds number 
swimmer model [5|, |20J . To incorporate catalytic activity, 
we take two enzymes, h and t, and place them at the head 
and the tail of the three-sphere device, respectively. The 
enzymes catalyze the following chemical reactions 



+ h^hF^ hQ- +G+ ^h + Q- +G+ 


, (1) 


- + P+ + t ^ tR- + P+ ^ t + S, 


(2) 



which would otherwise occur very slowly in the bulk. 
Since the head can be found in three distinct states of 
h, hF, and hQ~ , and the tail has two distinct states 
of t and tR~ , the combined system can have 6 distinct 
chemical states. The reactions ([Ij and ([2]) then multiply 
within these states and produce the chemical cycle shown 
in Fig. [H 

The second ingredient in the design is to take advan- 
tage of the presence of ionic components in the reac- 
tions and use electrostatic interactions to induce con- 
formational changes in the swimmer. If we make the 
middle sphere permanently negatively charged (see Fig. 
[2]), the temporarily charged states of hQ~ and tR~ will 
introduce electrostatic repulsion that can elongate the 
right and left arms, respectively. This allows us to map 
the 6 chemical states shown in Fig. [T] onto the 6 dis- 
tinct mechanochemical states shown in Fig. [21 These 
mechanochemical states can be represented in a 3D con- 
figuration space whose axes correspond to deformation of 
the right arm Ur, deformation of the left arm Ul, and 
complexation of h with F. The reaction routes in the 
chemical cycle in Fig. [1] can now be used to identify the 
corresponding transitions between the mechanochemical 
states in the configuration space, as shown in Fig. [2] 
The relevant transition rates are defined on the configu- 
ration space diagram in Fig. [2l I have chosen a notation 
where the rates of similar transitions are differentiated 
by a prime, with the general rule that k'^ > ka- These 
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Swimming Direction (J >0) 



FIG. 2. The 6 distinct mechanochemical states of the three- 
sphere swimmer in the configuration space. The conforma- 
tional states of the swimmer are determined by the projec- 
tion onto the space of the deformation of the left (Ul) and 
right (Ur) arms, while binding of h with F produces distinct 
chemical states that have identical conformations. All possi- 
ble transitions described in the chemical cycle in Fig. [l] are 
shown by arrows with attributed rate constants. 



choices are made depending on whether the presence of 
additional electrostatic charges hampers or facilitates a 
certain reaction 



211. 



We do not make assumptions about any specific or- 
dering of the various events involved in the cycle (see 
Figs. [T]and[l]). For example, state (1) can proceed to 
state (2) upon binding of R~ to t, and can also proceed 
to (3) if F binds to h. Any of these two can happen 
at a time in a stochastic manner, and only the average 
rates of these transitions are controlled. Moreover, not 
all of the transitions contribute constructively towards 
net swimming. For example, state (2) can proceed to 
state (4) upon binding of F to h, which will take the sys- 
tem further down towards completing a non-reciprocal 
mechanical cycle, but could also go back to state (1) if 
P+ binds to tR^ , which undoes the progress made by 
the previous transition towards net swimming. 

The dynamics of the device can be studied using a 
stochastic description of the conformational transitions 
and the relative population of the different states at sta- 
tionary state [20|. In stationary situation, each of the 
six mechanochemical states shown in Fig. [2] is populated 
with a probability pi (i = 1, • • • , 6), subject to the nor- 
malization condition pi + P2 + P3 + Pi + P5 + Pe = 1- 
The loops in the configuration space contribute prob- 
ability currents denoted as J, Ji, J2, and J3 in Fig. 
m Conservation of probability current leads to the fol- 
lowing nine equations between the probabilities and the 
currents: J + J3 = k^pe, J3 + Ji = kf[F]pi, J3 = 
^4P3, J2 - J ^ k'^Pb, J2 - Ji = kf[F]p2, J2 = kipi, 
Ji = k',[R~]p3 - kp[P+]p^, J = fc;[P+]p5 - kr[R-]pe, 
J — Ji = k'^[R-]pi — kp[P~^]p2- The above ten linear 
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FIG. 3. Dimensionless current j — J/k^ as a function of / = kf[F]/kT:,, r — kr[R~]/k5, and p — fcp[P"'"]/fc5, corresponding to 
a = 1.5, /3 = 1.5, 5 = 1.1, and h = \. (a) Michaelis-Menten behavior, for 7 = 1.1. (b) Sign change, for / = 1 and 7 = 2.7. (c) 
Nonmonotonic concentration dependence, for 7 = 1.1 and / -^ 00. 



equations can be solved for the ten unknowns (six prob- 
abilities and four currents) as functions of the concentra- 
tions {F\, {R-\, and [P+]. 

The six states in the mechanochemical configuration 
space of the swimmer correspond to only four indepen- 
dent conformations (Fig. [2]). To study the coupling be- 
tween the kinetics of the mechanochemical cycle and the 
hydrodynamics of the ambient fluid, we need to consider 
only the transitions between the distinct conformational 
states. Assuming that the extensions of the arms are 
small compared to their equilibrium lengths, it has been 
shown that the net swimming velocity of the device is 
controlled by the surface area A enclosed by the confor- 
mational states that make the (1) ?^ (2) <— (5) ^ (6) ^■ 
(1) loop (in Fig. [2]) and the net rate of sweeping this area 
J [20| . If we choose to make all the rates dimensionless 
by /c5, and in particular define the dimensionless current 
i = J /k^^ we can write the average swimming velocity of 
the device as 



gDAk^ . 



kf[F] kr[R- 



MP+]\ 






(3) 



where g is a geometrical numerical prefactor of order 
unity [22:], D is the diameter of the spheres, and L is 
the average total length of the swimmer. 

For simplicity, let us define the dimensionless variables 
a — k'^/kr, /3 — k'p/kp, 7 = ^4/^4, S — kyk^, and 
b — k^/k^, as well as the dimensionless concentrations 
/ = kf[F]/k5, r = kr[R-]/k5, and p = kp[P+]/k5. As 
argued above, we expect a, f3, 7, and S to be numeri- 
cal factors (not too much) larger than unity. Solving the 
above system of equations, we find that detailed-balance 
holds (i.e. j = 0), when 



7(a/3 -6) + (a/3 -"fS){f +p + ar) = 0. 



(4) 



By tuning the concentrations /, r, and p away from val- 
ues that satisfy the above equality, we can drive the sys- 
tem away from equilibrium and obtain a nonvanishing 
swimming velocity {j ^ 0). 



The concentration dependence of j reveals a number 
of interesting features in the device, as shown in Fig. 
[3l Figure [3^ shows that the dependence of j and conse- 
quently the swimming velocity on the concentration of F 
adopts a Michaelis-Menten form. This behavior, which is 
also observed in molecular motors [15|, is inherited from 
the (multi-step enzymatic) reaction kinetics. Figure |3Jd 
shows that for sufficiently low concentrations of F and 
relatively large values of 7, it is possible to have a situa- 
tion where j changes sign, which leads to the reversal of 
the swimming velocity. This can be understood from the 
condition of detailed-balance in Eq. ((4]) as j is propor- 
tional to its left hand side. Figure [3fc shows that when 
the concentration of F is at its saturation limit (see Fig. 
[3li) the dependence of j on the concentrations of R^ and 
P"*" is nonmonotonic, and suggests a strategy to optimize 
the swimming velocity. 

There are many efficient enzymatic reactions that 
could be used in practice in such a setup. For example, 
Eq. Ill can represent the conversion of carbon dioxide 
{F — CO2) (plus water) into bicarbonate {Q~ = HCO;^) 
and proton (G"^ = H+) catalyzed by carbonic anhydrase 
{h — CA) [23|. An example for Eq. ^ could be su- 
peroxide anion (i?" = 02-~) being scavenged by su- 
peroxide dismutase {t = SOD) by consuming a proton 
(P+ = 11+) and producin g hy drogen peroxide and oxy- 
gen {S = ^HzOa -I- iOa) [ll]. Looking at Figs. [T]and 



we find that in this example k^ = kcA 
[231, while kr = fenn ~ 3 x 10^ 



6 X 10^ s"^ 
23. While 



'SSOD ~ J X lU" M~ s 

we need to know all of the rates to be able to perform a 
full calculation, it is possible to make a rough estimate 
by looking at typical optimal values of j '^ 10 ~^ from 
Fig. [21 and the above estimate for k^. Using Eq. ([3|) 
with D ^ 1 nm, A '^ 1 nm^, and L '^ 10 nm, we find an 
estimate of ^ ~ 100 nm s""'^, which is comparable to the 
efficient molecular motors that use protofilament tracks 
for their locomotion. In this molecular realization, these 
enzymes could presumably play the role of the spheres 
themselves, and could be attached to one another using 



unstructured tails or covalently bonded polymeric teth- 
ers. It might also be feasible to attempt a larger scale 
setup with colloidal spheres covered with adsorbed en- 
zymes and connected with polymeric tethers. 

In practice it might be difficult to fully suppress other 
transitions, and we might have to deal with more compli- 
cated configurational spaces (e.g. involving other tran- 
sition routes between these states) as well as the back 
reactions. Moreover, realistic chemical reactions might 
have more complicated kinetics (such as the four-stage 
kinetics of the reaction of superoxide anion by SOD [23[ 
that we have simplified into two stages in the above ex- 
ample) . Such complications can be taken into account in 
the calculation of the average swimming velocity using a 
similar strategy [25|, although they are not expected to 
change the qualitative behavior of the system provided 
the detailed-balance is broken. 

In addition to affecting the conformational transitions, 
thermal fluctuations will also tend to randomize the ori- 
entation of the swimmer, on time scales larger than the 
rotational diffusion time t^ ~ rjL'^ j (k^T) , where rj is the 
viscosity of water, fee is the Boltzmann constant and T is 
the temperature. In addition to the rotational diffusion 
of the swimmer, fluctuations in the density of the fuel 
molecule could also induce anomalous fluctuations in the 
magnitude and direction of the velocity. The combina- 
tion of these effects could lead to a variety of different 
behaviors, similar to those found in the context of self- 
phoretic swimmers [26]. We note that a related enhanced 
diffusion of active enzymes has recently been observed ex- 
perimentally 27| . The dependence of the velocity of the 
swimmer on concentration of the fuel in Eq. ([3]) also sug- 
gests that such artificial swimmers could be guided using 



chemotaxis 18|, |28|, |2 



In conclusion, I have proposed a design for a molecu- 
lar scale low Reynolds number swimmer that can use the 
free energy gain and products of a nonequilibrium chem- 
ical reaction and undergo conformational changes that 
are non-reciprocal and therefore can lead to net swim- 
ming under low Reynolds number conditions. While the 
present analysis has used a minimal kinetic model that 
can achieve swimming at such small scales, the scheme 
can be readily generalized to help analyze any realistic 
mechanochemical cycle corresponding to more complex 
molecular devices that could be synthesized. 
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